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A simple first-principles method is proposed to describe propagation of tightly bound excitons 
via its kinetic kernel. By viewing the exciton as a composite object (internally consisting of local 
particle and hole in Wannier orbitals), the corresponding kinetic kernel encapsulates the detailed 
information of the propagation and decay of the exciton, applicable for all binding energies. Case 
study of LiF gives three branches of exciton dispersion, which we verified via inelastic X-ray scat- 
tering. The proposed real-space picture not only offers computationally inexpensive access to the 
full exciton dynamics, suitable even in the presence of surfaces and impurity scattering, but also 
provides intuitive physical understanding necessary for the optimization of exciton propagation in 
semiconductors, their nanostructures, and strongly correlated oxides. 



PACS numbers: 71.15.Qe, 71.35.-y, 78.20.Bh, 78.47.-p 

Exciton, a low-energy charge excitation involving a 
bound particle-hole pair, is one of the most utilized ex- 
citation of materials in practical applications. It plays 
important roles in many semiconductor applications, for 
example, in photovoltaic and other PN-junctions related 
devices [THS]. One essential issue in real applications is 
the efficiency of the propagation of excitons and the asso- 
ciated decay / recombination processes [THS]. Currently, 
excitons are typically studied via perturbation theory 
(so-called Bethe-Salpeter equation) [7HTT], or via spe- 
cially tuned approximation within the time-dependent 
density functional theory |llH14j . Particularly, nonim- 
mune to the high computational expense, the emphasis 
has been only on the excitonic spectrum in the optical 
limit, most relevant to the creation of zero-momentum 
exciton, leaving the essential issues of exciton propaga- 
tion mostly unaddressed. 

The recent advances in X-ray scattering that allows 
mapping out the full momentum-energy dispersion of the 
excitons |15j further highlights the need to have a compu- 
tationally inexpensive method, not only to describe the 
exciton propagation, but also to give an intuitive physical 
picture and a deeper understanding. This is particularly 
so given that more and more of the potential functional 
materials make use of the strongly correlated nature of 
the electrons, which require treatment of the many-body 
problem beyond the scope of perturbation theory. It is 
thus important and timely to develop a first-principles 
theory targeting specifically the exciton propagation. 

In this letter, we propose a general approach that fo- 
cuses on the essence of the kinetics of exciton. By inte- 
grating out the higher-energy virtual pair-breaking pro- 



cesses, the exciton can be treated as a single composite 
object, instead of a particle-hole pair, and its kinetics 
can then be formulated in real-space via an excitonic ki- 
netic kernel that serves the role similar to the hopping 
integral. In addition, the kinetic kernel encapsulates the 
information of exciton decay process and the necessary 
energy dependence to react to all binding energies. For 
systems where excitonic effect is large, namely the bind- 
ing of exciton is strong, this real-space description, to- 
gether with the reduction of the problem from two-body 
to one-body, is significantly more efficient for computa- 
tion (easily by several orders of magnitude). In addi- 
tion, the separation of non-local kinetic from the local 
interacting problem allows an accurate treatment of the 
local many-body problem. Taking the well known LiF 
as a case study, our results give three branches of exci- 
tons, which we verified experimentally via inelastic X-ray 
scattering. The proposed approach offers a simple and in- 
tuitive physical picture easily applicable to systems with 
surfaces (or other boundary conditions), impurity scat- 
tering, or strong many-body interactions. Therefore, it 
should facilitate greatly understanding of exciton prop- 
agation in all semiconductors, their nanostructures, and 
numerous strongly correlated functional materials. 

The general idea of our approach is to use local exci- 
tons as real space basis and encapsulate their propagation 
into an excitonic kinetic kernel, similar to the hopping 
integral. This real-space approach makes sense because 
binding of the particle and hole is only significant in short 
range and thus the size of the exciton with practically rel- 
evant binding should be reasonably small. To this end, 
one needs to 1) define a "local" region, in which the local 



excitons would be treated accurately enough, and 2) con- 
struct the kinetic kernel of the local excitons which gives 
a physical understanding of the exciton propagation and 
allows an inexpensive computation of the exciton spec- 
trum in the full momentum / energy phase space. As a 
bonus, the separation of the kinetic processes from the 
local interacting problem would also allow a more accu- 
rate treatment of the local many-body physics. The first 
step is strongly case dependent, in both the physical size 
of the local region and the required proper treatment of 
local interaction (e.g.: local multiplets in strongly cor- 
related materials) . The case study of LiF below gives a 
detailed example of this step. The second step, on the 
other hand, is quite generic and will be the focus of this 
study. 

Let's define the excitonic kinetic kernel T through the 
one-body equation of motion of the exciton (in matrix 
notation) 



L[H] = L[H L ] + L[H L ]TL[H], 



(1) 



given the particle-hole correlation function L correspond- 
ing to the local Hamiltonian Hl, and the full Hamilto- 
nian H with additional coupling beyond the local region. 
Note that the excitonic degree of freedom b RN , as a ba- 
sis for I/ra^r/ap and Tr^r/a^, is defined in the direct 
product space of the particle p and hole h orbitals at 
the same local region indexed by R: b RN = c^cra, 
since only strongly bound excitons smaller than the local 
region are of interest. Eq.[l] makes the physical interpre- 
tation of T obvious: T is the kernel responsible for the 
propagation (present in L[H]) of local excitations (given 
hyL[H L }). 

This rigorous definition becomes practically useful 
with the simplification: 

T =L~ 1 [Hz] — 

={L^[G[H L ]] I[H L ]) (L^[G[H}\ I[H\) (2) 

^L^[G[H L \\-LZ l [G[H\\- 

Here L (HN, H'N'; t, t') = G(Rp,R'p';t,t') * 
G(R'h',Kh;t',t) is the "bare bubble" without in- 
teraction (vertex) at the two-particle level (even though 
it does contain the fully interacting G[H}.) This sim- 
plification of the exciton self-energy I[H] by its local 
counterpart I[Hl\ is similar to that previously employed 
with the dynamical mean field approximation [16 . It is 
physically reasonable because of the general short-range 
nature of the self-energy. Particularly, at the lowest 
order approximation of /, employed in practically all 
first-principles solution of the Bethe-Salpeter equation, 
I[H] = I[Hl] and this simplification is exact. Once T is 
obtained from Eq. [2] the full spectrum of L[H] can be 
obtained effortlessly from Eq. [T] with several orders less 
computation effort than the two-body Bethe-Salpeter 
equation (BSE). 
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FIG. 1: (color online). Three ways to define local regions 
via unit cells respecting to (a) translational symmetry, (b) ro- 
tational symmetry, and (c) both translational and rotational 
symmetries. Following (c), an effective F-s' Wannier orbital 
(d) composed of symmetric superposition of Li 2s orbitals is 
compared to a F 2p Wannier orbital (e). (f) Shape of the 
exciton from product of Li 2s and F 2p gives the form factor 
(g) in momentum space. Boxes denote the conventional unit 
cell in (d), (e), and (f), and the corresponding 1st Brillouin 
zone in (g). 



We now use the prototypical LiF to demonstrate the 
simple method and the deep insights it can offer. Crys- 
tal LiF has a simple NaCl structure with ionic electronic 
structure, with a 14.2 eV band gap p~7] defined by oc- 
cupied F p bands and a unoccupied Li s band. The 
physics of excitons in LiF is thus dominated by these 
low-energy bands. Furthermore, the strong ionicity of 
the system implies a dominating short-range attraction 
between a hole in the F p orbitals and a particle in the 
nearest neighboring Li s orbital. We will therefore define 
a local region containing only one F for simplicity. The 
straightforward choice is a small unit cell containing one 
F and one Li, as illustrated in Fig. [ija). This choice, 
however, does not build in the rotational symmetry in 
the local region. Fig. |TJb) gives an alternative choice of- 
ten employed in the cluster model [TS]. While including 
the rotational symmetry, the local Hilbert space is, how- 
ever, overcomplete in that one F is paired with four Li 
(actually six Li in 3D), and thus cannot be used to build 
a crystal under translational symmetry. Fig. [ljc) gives 
the best choice of the local region, which contains one F 
and effectively one Li from combining all its neighboring 
sites. This choice respects both the translational and ro- 
tational symmetries of the system and is thus ideal for 
the physical description. 

In practice, the low-energy Hilbert space of such a local 
region can be constructed by making use of the gauge 
freedom of the Wannier function [HO HO] • Fig. [ljd) shows 
the resulting Li s Wannier orbital centered at the F site, 
as a symmetric superposition of atomic s orbitals. Recall 
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a similar construction [5T] in the well-known Zhang-Rice 
singlet in the cuprate. Effectively, the local region is just 
a "super F atom" containing three p orbital (cf: Fig. We)) 
and one s' orbital [TS] in a simple fee lattice, and the 
tightly bound excitons can now be viewed effectively as 
intra- "atomic" Frenkel excitons. The excitonic Hilbert 
space N is then simply spanned by three local excitons 
made of a particle in s' orbital and a hole in p x , p y , or 
p z orbital. Fig. [ljf) also demonstrates the probability 
amplitude of the exciton in real space, obtained from the 
product of the particle and hole orbitals 22J. 

With the local Hilbert space defined, excitonic kinetic 
kernel can be easily computed via Eq. [2j As a sim- 
ple illustration, the fully dressed Green's function is ob- 
tained approximately from self-consistent DFT Hamilto- 
nian with a +5 eV shift in the on-site energy of s' orbital 
to counter the self-interaction error. The corresponding 
band structure contains a 14.2 eV band gap, consistent 
with the GW calculation [T71 12"5] . We then compute the 
corresponding Lq[Hl\ and Lq[H] from the standard con- 
volution [53] in Wannier function basis. The resulting 
Lq[Hl] contains trivially a pole at ujq = lu p — u>h = 21.9 
eV, corresponding to the energy difference between the 
particle and the hole. (Since all three "bare" excitons 
are symmetry related, Lq[Hl] is diagonal and N inde- 
pendent.) Turning on the local binding interaction of 
Ub ~ J dxdx'V S p(x,x')n s (x)np(x') ~ 7 eV simply lowers 
the pole to uj exciton = u - U B — 14.9 eV in L[H L ] (ig- 
noring the fine multiplet splitting.) Note that our Ub is 
with respect to the binding energy of local bare exciton 
and is thus much larger than the binding energy between 
Bloch-like orbitals, typically referred in the context of 
Wannier exciton. 

As shown in Fig. [2j the resulting excitonic kinetic ker- 
nel obtained from Eq. [2] contains several interesting char- 
acteristic features. First, the kinetic kernel is strongly 
energy dependent. This is physically necessary because 
T needs to generate the correct [/^-dependent kinetics of 
excitons without the knowledge of Ub to be introduced 
in L[Hl]. Particularly, T hosts a l/(w — ioo) = ^/Ub 
decay in all its components. Since T is like the hopping 
integral, this indicates a heavier mass (reduced hopping) 
for excitons with stronger binding (lower energy) . This is 
the exact behavior in the strong coupling limit, in which 
the exciton can only propagate by virtually moving the 
particle first (paying the binding energy U B in the inter- 
mediate state,) followed by the hole (to reform the pair), 
a 2nd order process giving T ~ 2t p th/UB oc 1/Ub- (Here 
t p /th refers to the hopping of the particle/hole.) 

Second, the kinetic kernel has imaginary part spanning 
over a wide energy range, representing excitonic decaying 
channels. In fact, the energy span of the imaginary part is 
exactly that of the Landau continuum from unbound ex- 
citations. This is not surprising because plugging L[H L ] 
with unbound excitons (Ub = and L[Hl] = Lq[Hl\) 
into Eq. [I] must, by construction, recover exactly the 
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FIG. 2: (color online). The real (solid lines) and imag- 
inary (dashed lines) parts of two-particle kinetic kernel 
T(RN, R'iV) for (a) R' — R=(0,0.5,0.5), (b) R' — R=(0,0,1), 
and (c) R' — R=(0,0,0) for N = p x , p y , and p z . Virtual pro- 
cesses integrated in off-site (d) and on-site (e) components of 
T. 



TABLE I: ReT(RiV, R'iV'; u exciton = 13.4 eV) in meV. 
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Landau continuum in the resulting L[H] = Lq[H], and 
the information of the continuum can only come from T. 
(L [Hi,] has no knowledge of the dispersion.) Physically, 
this reflects the fact that the definition of Trjv,r/ N' y i a 
Eq. [T] integrates out all the virtual processes with un- 
bound particle and hole traveling independently from R' 
to R, as demonstrated in Fig.[2](d). Therefore, the T de- 
scribes the exciton kinetics exactly not only in the strong 
binding limit, but also in the weak binding limit. 

Third, the kinetic kernel has a significant negative di- 
agonal elements [cf. Fig. |2fc)] that reduces the local 
exciton energy to ui exc iton ~13.4 eV (estimated from 
uj — LOexdton — ReToiv,oJV (w) = 0.) This renormalized 
exciton energy represents the average energy of the prop- 
agating exciton in the crystal. It is important to distin- 
guished this renormalization from the binding of local 
exciton, as it originates from virtual kinetic processes il- 
lustrated in Fig. ^e). For large Ub, the lowest order of 
such processes would give T ~ — (tp + tfy/Us for each of 
the 12 possible intermediate de-paired states in the FCC 
lattice. 

Before computing the full L from Eq. [T] let's first ex- 
amine the physical insights Re T offers, by considering it 
as an effective hopping integral (cf. Table [I]). Regarding 
and Re TR iR '^ R (aj exc i ton ) as diagonal and off- 
diagonal elements, a simple diagonalization gives an esti- 
mated exciton dispersion containing three exciton bands 
shown in Fig. [3]ja). Notice a striking similarity of the 
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FIG. 3: (color online), (a) Exciton dispersion from diago- 
nalizing HeT(cu = 13.4 eV) (solid curves), compared to £(q) 
colored by red. (b) Band structure of Wannier orbital for 
Li s (E > 0) and F p (E < 0). (c) Selected k paths on 
the xy plane, denoted by <j>. (d) Exciton dispersion along se- 
lected paths, with the longitudinal component colored in red. 
Dashed lines denotes the zone boundaries in (c). 



FIG. 4: (color online), (a) Imaginary part of dynamical 
linear response function % m selected directions (c.f.:Fig. |3j, 
together with exciton dispersion from ReT at ui — 13.4 eV 
(solid curves). The theoretical intensity for q < 1.5 is reduced 
by 2.6 to mimic screening by the plasmon pole, (b) Exciton 
dispersion measured by inelastic X-ray scattering in the same 
directions. 



dispersion to that of the F-p bands in Fig. |3jb), beside 
a smaller bandwidth. This is understandable since the 
excitons here have the same symmetry properties as the 
F-p orbitals [c.f. Fig. |T](e) and (f)], thanks to the zero 
angular momentum of the F-s' orbital. Consequently, 
T ~ 2(t p /V r s)th simply resembles th with a 2t p /UB re- 
duction that wipes out all longer range hoppings. (The 
comparison appears up-side-down because t p < 0.) Ob- 
viously, the standard effective mass description [25] US] 
of Wannier exciton propagation is seriously inadequate 
here. In essence, the exciton binding always enhances 
the tight-binding nature of the excitons, making them 
propagate more like Frenkel excitons. 

Next, the complete description of exciton propagation 
and decay are computed from inserting T into Eq. [T] In- 
deed the resulting dispersion is quite close to the above 
simple estimation except: 1) Near the bottom of the exci- 
ton bands around F point, the actual dispersion is weaker 
than the estimation, due to the above mentioned decay 
of T(lj) at lower energy. 2) Near the top of the exci- 
ton bands, the exciton peak is physically broadened, re- 
flecting a shorter lifetime from decaying into the Landau 
continuum encapsulated in ImT. 

Our results of exciton propagation leads to an excit- 
ing realization that all three branches of the excitons 
can be observed via inelastic X-ray scattering. Fig. J3jd) 
shows the exciton bands along different paths in the Bril- 
louin zone, weighted by their longitudinal (|| q) contribu- 
tion. Typically, one would expect that out of the three 



branches of the exciton, only the longitudinal mode can 
be observed experimentally, given the X-ray matrix ele- 
ment M$ = /dxexp" iq ' x (x|0j3)(0ft,|x) [22 . However, 
because of the zone folding, hybridization with trans- 
verse modes becomes significant near the zone boundaries 
in most off-high-symmetry directions. Consequently, the 
transverse modes receive significant spectral weight out- 
side the first zone, allowing direct experimental observa- 
tion of all three branches of excitons. 

We verify the theoretical results by performing inelas- 
tic X-ray scattering along different directions between 
(100) and (110) [24]. Indeed, the measured spectra show 
good agreement with the theoretical response function, 
X(q, w ) = EjviV' M N(<l)L NN >(q)M^,(q), despite the 
resolution-related broadening in the experimental data. 
(A spectral weight reduction for q < 1.5 due to plasmon 
screening is included in L via additional Hartree contri- 
bution [27 . This improves the agreement in the spectral 
weight, but has little effects on the exciton dispersion 
of interest here.) As soon as q deviates from the high- 
symmetry direction, a clear "break" in the dispersion is 
observed, for example at q ~ 3.5(27r/a) for 9 — 7.5° and 
at q ~ 4(27r/a) for 37.5°, which actually corresponding to 
switching from one exciton band to another, as indicated 
in the theoretical results. 

Finally, several important aspects of our method is 
worth mentioning. First, the separation of the non-local 
kinetics from the local binding of exciton in our frame- 
work allows a more accurate many-body treatment of the 
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local excitation in L[Hl] (say via exact diagonalization 
or quantum Monte Carlo) before solving Eq. [I] This is 
particularly important in strongly correlated systems, in 
which local interaction is strong enough to induce sig- 
nificant multiplet splitting [57]. Furthermore, such sep- 
aration of the problem significantly reduces the compu- 
tational expense of exciton dynamics, as the proposed 
evaluation of excitonic kinetic kernel T is purely non- 
interacting, while the potentially expensive interacting 
problem is only local. On top of that, our use of Wannier 
function as basis, and the implicit projection of non-local 
particle-hole pair from our formulation further reduce the 
computation greatly. We estimate that our method is at 
least several orders of magnitude cheaper than the stan- 
dard Bethe-Salpeter equation where a ~ 10 8 matrix el- 
ements of the electron-hole interaction [53] are needed 
in comparison with ~ 10 4 matrix elements. A similar 
reduction (by ~2 orders of magnitude) is achieved in 
the matrix elements of the form factor of the experimen- 
tal probe. Finally, the real-space formulation employed 
in this study is particularly suitable in treating exciton 
scattering from surfaces, impurities, or other boundary 
conditions. All together, these significant advances make 
this method ideal for future studies of exciton propaga- 
tion and decay in modern functional materials, like doped 
semiconductors and strongly correlated transition metal 
oxides, as well as their nano-structures, thin films, and 
heterogeneous interfaces. 

In summary, we propose a new real-space theoreti- 
cal framework to describe the propagation and decay of 
tightly bound excitons. By explicitly integrating out the 
intermediate de-pairing processes, a simple one-particle 
description of exciton is given by the effective exciton 
kinetic kernel, which gives intuitive physical picture of 
both mass enhancement and decay into unpaired particle 
and hole. The proposed method is exact in both weak 
and strong binding limit. Furthermore, the separation 
of non-local kinetic from the local binding of excitons 
allows a more accurate treatment of local many-body 
problem beyond perturbation theory, particularly nec- 
essary for strongly correlated materials. Together with 
the Wannier function approach, the method is compu- 
tationally much more affordable than the standard two- 
particle Bethe-Salpeter equation method, by easily sev- 
eral orders of magnitude. In the prototypical LiF, our 
theory predicted three exciton bands that can be mea- 
sured experimentally over a large momentum. We verify 
this via inelastic X-ray scattering and found good agree- 
ment. Our real-space approach can incorporate straight- 
forwardly scattering via impurities or boundary condi- 
tions, and thus will find broad applications in the study of 
nano-structures, thin films or surfaces of semi-conductors 
and strongly correlated Mott insulators. 
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